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ABSTRACT 

We present a detailed analysis of the iron Ka line variabilities in NGC 
4151 by using long ASCA observation data obtained in May 1995. Despite the 
relatively small amplitude variations in the continuum flux, the iron Ka line flux 
and profile show dramatic variations. Particularly the line profile changes from 
single peak to seeming double peaks and back in time scales of a few 10 4 sec. 
The seemingly double-peaked profiles can be well interpreted as line emission 
from a Keplerian ring around a massive black hole. An absorption line at around 
5.9 keV is also marginnaly detected. We discussed current Fe K line models, 
but none of them can well explain the observed line and continuum variations. 

Subject headings: black hole physics — galaxies: active — galaxies: individual 
(NGC 4151) — line: profiles — X-rays: galaxies 
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1. Introduction 

Broad Fe K fluorescence line profiles, detected in many AGNs (Tanaka et al. 1995, 
Nandra et al. 1997, Wang et al. 1999a, Nandra et al. 1999, and references therein), are 
thought to arise from the accretion disk near the center massive black hole of AGNs. The 
line profile, thus, carries important information about the distribution and kinematics of 
optically material as well as about the space properties in the vicinity of the postulated 
super massive black hole. Doppler and gravitational shifts would imprint characteristic 
signatures on the line profile which map the geometric and dynamical distributions of 
matter surround the black hole. 

Additional information concerning the geometry of X-ray source and matter in the 
active nucleus could in principle be derived by studying the rapid variability (on time 
scales of several 10 4 s or less) of the line profile, intensity and their relationship with the 
continuum variations. So far, rapid variability in Iron K line has been detected in NGC 
7314 (Yaqoob et al. 1996), MCG -6-30-15 (Iwasawa et al. 1996, 1999), NGC 4051 (Wang 
et al. 1999b) and NGC 3516 (Nandra et al. 1999). 

The observed Fe Ka line rapid variabilities in different targets or during different 
observations of the same target are quite different. Yaqoob el al. (1996) presented the 
evidence for rapid variability of the Fe K line profile in the narrow-line Seyfert galaxy NGC 
7314, which is consistent with a disk-line of constant equivalent width superposed on a 
constant flux narrow line (presumably from the torus). Iwasawa et al. (1996) discovered 
that, during the ASCA observation on MCG -6-30-15 in 1994, when the source is bright, 
the Fe K line is weak and dominated by the narrow core, whilst during a deep minimum, 
a huge red tail appears. The intensity of broad Fe K line correlates inversely with the 
continuum flux on the time scales of several 10 4 s. And Iwasawa et al. (1999) discovered 
that during the ASCA observation on the same target in 1997, the Fe Ka line of a flare 
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phase peaks around 5keV and most of its emission is shifted to below 6 keV with no 
component detected at 6.4keV. While for NGC 4051, during the ASCA observation in 1994, 
the equivalent width and the width of Fe K line correlate positively with the continuum 
flux (Wang et al. 1999b), which shows an opposite trend with the ASCA observation on 
MCG -6-30-15 in 1994. The rapid variabilities of Fe Ka line profiles in NGC 3516 are also 
quite irregular (Nandra et al. 1999, Wang et al. 2000). 

The nearby Seyfert 1.5 galaxy (z = 0.0033), NGC 4151, was first established as a 
bright X-ray source some thirty years ago (Gursky et al. 1971). Since that time, this 
archetypal low-luminosity active galaxy has invariably been considered as a prime target by 
all major X-ray astronomy missions. In 1995, NGC 4151 was observed by ASCA from May 
10 to 12 (Leighly et al. 1997). The net exposure time of the observation is about 200 ks. 
Because of the enhanced sensitivity provided by ASCA and the much longer exposure time, 
the obtained Fe Ka line profile has very high data quality (Wang et al. 1999a, here after 
W99). In this letter we report the discovery of unique and rapid variabilities of the iron line 
profiles in NGC 4151 during the ASCA observation in 1995. 



2. Fitting the X-ray continuum 

Only the SIS data are used in this paper to study the Fe Ka line variabilities. Detailed 
description of the ASCA observation and data reduction can be found in W99. Figure 
la shows the background subtracted light curve (0.4 - 10.0 keV) with 256 s binning. 
Significant variability is absent on time scales of less than 1000s, but moderate variations 
can be seen on time scales of several 10 4 s, exhibiting a continuous slow drifting behavior 
(e.g. Lawrence 1980). The light curve shows an intensity dip at 60 ks from the beginning, 
which lasts about 50 ks, and two neighboring flares at about 150 ks, for which each lasts 
about 30 ks. The intensity ratios of the flares to the dip are around 0.8. Apart from these, 
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the count rate also increases slightly during the whole observation. 

The time-averaged Fe Ka line profile from this observation has been extracted by W99. 
The data quality of the line profile is one of the best that could be found in literatures 
(see Figure la in W99). In order to analysis the rapid variability properties of the Fe Ka 
line profiles, we divide the whole observation into five time intervals (see Figure 1, 1-1 to 
1-5). Following Weaver et al. (1994) and W99, we fit the underlying continuum of the each 
segment in the 1.0 - 4.0 and 8.0 - 10.0 keV bands (to exclude the broad iron line region) 
with a model which consists of a dual absorbed power law plus some fraction (the best-fit 
value for the time averaged spectrum is ~ 5%) of the direct continuum scattered into our 
line of sight and absorbed only by the Galactic column of 2 x 10 20 cm -2 . We do not include 
a Compton reflection component in the spectral fits because of its small impact on the 
ASCA spectra. During the fitting, we fix the power law index and the flux of the scattered 
component at their best-fit values obtained from the time-averaged spectrum. Same to the 
time-averaged spectrum, the spectra of these five segments can also be well fitted by this 
model. We show the results in Table 1. 

The power law indices for individual segments are all well consistent with the average 
value, i.e, there is no significant change in the continuum spectrum. However, the absorption 
columns density and the covering factor of the dual absorber showed significant variations 
(table 1), not correlated with variations of the continuum flux. 

3. The Fe Ka line profiles 

Figure 2 shows the Fe Ka line profiles for the whole observation and the five time 
intervals. The profiles are obtained from the ratio, the data divided by the best-fit 
continuum (folded through the detector response), multiplied by the power-law (original 
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form), thus are detector effective area corrected, and also independent from the model used 
to fitting the lines. The time-averaged line shows a strong narrow peak around 6.4 keV and 
a huge red wing extending to ~ 4.5 keV. Apart from a similar narrow core in the 6.4 keV, 
the line profile apparently differs by one interval to the other. The line profile of 1-1 also 
shows a huge red wing, extending to about 4.0 keV with a possible red peak at around 4.9 
keV; while for 1-2, the red wing is much weaker. The line profile of 1-5 is quite similar to 
the time averaged one. For 1-3, besides the strong peak at 6.4 keV, there is also a seeming 
peak at around 5.4 keV. We noticed that between the two peaks there seems a valley in 
between which indicates that the line is double-peaked. The line profile of 1-4 is similar to 
that of 1-3, seemingly double-peaked, but much weaker. We conclude here that the Fe Ka 
line profiles of NGC 4151 show clear and obvious variations during the ASCA observation. 

Initially the line profiles are fitted with two gaussians and the results is given in Table 
2. A narrow gaussian at around 6.4 keV and additional broad component can well fit the 
line profiles, and the results clearly present the variations of the line profiles and equivalent 
widths. 

We also fit the Fe K line profiles with a model consist of a narrow component and a 
disk line (Fabian et al. 1989). The narrow core is likely emitted by material far away from 
the center black hole, such as torus, thus we would not expect substantial change in the line 
flux within the duration of the observation. During the fitting, we fix the flux of the narrow 
core for the five time intervals at the best-fit value from the time-averaged spectrum (which 
has a equivalent width of 70 eV according to the time-averaged spectrum). We also fix the 
disk-line energy at 6.4 keV in the source rest frame (W99), and fix the inclination angles of 
the disk at the best-fit value from the time-averaged spectrum. In fact, when we free the 
inclinations for the five time intervals, the best-fit values are all well consistent with that of 
the time-averaged spectrum (27 degree). As W99 did, the corresponding Fe K/3 and Ni Ka 
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components (disk-line plus narrow core, see W99 for detail) are also included in the fit. The 
model can fit the line profiles quite well as the two Gaussians model do, and the results of 
the fit are shown in Table 3. We did not adopt the Model B in W99 because the scattered 
disk line component should be constant on such a short time scale, which is inconsistent 
with the rapid Fe K line variabilities we detected. 

Nandra et al. (1999) claimed that they discovered a red shifted iron resonant absorption 
line at around 5.9 keV in the Fe Ka profile of NGC 3516. We also add an additional 
absorption line component to the model, and find that including such a line can significantly 
improve the fit for 1-4 (5\ 2 = 7). However, the absorption line is not required for other 
segments (5x 2 < 3.0). 

4. Discussion 

4.1. The seemingly double-peaked line profiles 

Massive black holes are generally thought to exist at the centers of active galaxies (Rees 
1984), but unambiguous identification of a black hole has been impeded by lack of evidence 
for the strong-field relativistic effects expected in the vicinity of black holes. Because of 
its unique anticipated profiles, the Fe Ka fluorescence line has long been considered as an 
observational key to the study of the innermost regions of active galactic nuclei (AGN). 
Fabian et al. (1989) and Laor (1991) pointed out that the line profiles from the accretion 
disk should be broad, skewed and double-peaked (dependent on the inclination of the disk 
and the line-producing region in the disk), which are characteristics of the Doppler effects of 
the accretion disk and the strong gravitational field in the vicinity of the central black hole. 

The Fe Ka line profiles detected in many AGNs have been proved to be skewed and 
extremely broad (Tanaka et al. 1995, Nandra et al. 1997, Wang et al. 1999a, Nandra et al. 
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1999, and references therein), consistent with the disk line model, so give the ever strongest 
evidence for the presence of massive black holes and accretion disks in the center of AGNs. 

However, some alternative models, such as Comptonization and jet, can also explain 
the broad and skewed Fe Ka line profiles. Those models do not require strong gravitational 
field and so give no support to the presence of massive black holes (Fabian et al. 1995, 
Misra & Sutaria 1999). But different to the disk line model, those models predict a smooth 
line feature, but not a double-peaked one. Though some indirect evidences have shown 
that they are not satisfactory (Fabian et al. 1995, Misra 2000, Reynolds & Wilms 2000), as 
Misra & Kembhavi (1998) have pointed out, we infer that the detection of a double-peaked 
line profile is the key to confirm the disk line model for the observed Fe Ka lines, and so on 
confirm the detection of extremely strong gravitational field. 

In this paper we find that the Fe Ka line profile of NGC 4151 is seemingly double- 
peaked during phase 1-3 and 1-4 in Figl. Because an additional resonant absorption line at 
5.9 keV has also been detected in the line profile in 1-4, we only focus on the profile of 1-3 
in this section. We show the seemingly double-peaked line profile of 1-3 in detail in Fig. 3, 
and the best-fit two gaussians and disk-line models are also plotted. 

The seemingly double-peaked line profile strongly suggest that the line arises from 
materials with disk-like geometry. Additional strong gravitational field is also needed to 
explain the red shift of the profile. Others models for the broad line profiles, such as 
Comptonization or jet, can not be responsible for a double-peaked line profile, and disk-like 
outflows or inflows with high velocity can not be responsible for the red shift of the line 
profile. The best disk line fits show that the line mainly arises from a narrow ring between 
10.8 R g and 18.6 R g . 

Currently, perhaps due to the limited energy resolution and data quality, obvious 
double-peaked Fe Ka line profile has never been detected in any active galactic nuclei. 
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We noticed that the line profile of MCG -6-30-15 in Fig. lb of Wang et el. (1999a) looks 
double-peaked, but the plot was folded through the effective area of the detectors which 
rises sharply towards the low energy, so did not give the true profile of the line. While 
in plottings corrected for detectors' effective area, we can not see obvious double peaks 
(see Fig. 7 of Iwasawa et al. 1996 and Fig. 6 of Fabian et al. 2000). We hope that the 
new generation X-ray satellites, such as XMM, can show us clearly double-peaked iron line 
profiles of AGNs. 

4.2. The rapid variabilities of the Fe Ka line profiles 

Additional to the double-peaked line profiles, the rapid variabilities of the Fe K 
line profiles in NGC 4151 are also detected in this letter. We divided the whole ASCA 
observation of NGC 4151 in 1995 into five time intervals, and found that, despite the 
relatively small amplitude variations in the continuum flux, the Fe Ka line flux and profiles 
show dramatic variations: during the first time interval 1-1, the iron line profile is broad 
and skewed, showing a strong peak at 6.4 keV, and a huge red wing extending to ~ 4.5 
keV, similar to the time-averaged line profile; after several 10 4 s, the huge red wing becomes 
very weak (1-2); for 1-3, except the strong peak at 6.4 keV, there is also a seeming peak 
at around 5.4 keV; the line profile of 1-4 is similar to that of 1-3, but much weaker; and 
finally, the line profile change back to single peak (1-5). The line equivalent widths of the 
five intervals are 375, 148, 329, 168 and 300 eV respectively, showing obvious variations. 

The rapid variabilities of Fe K line profiles discovered previously have suggested that 
the X-ray continuum flux is generated by magnetic flares above the accretion disk, and at 
any given time there are only a few flares moving around above the disk (Iwasawa et al. 
1996, Fabian 1997, Wang et al. 1999b, Iwasawa et al. 1999, Wang et al. 2000). Flares with 
different location and different motion can give different fluorescence line profiles. 
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Can the magnetic flare model also explain the observed rapid variabilities of Fe Ka 
line in NGC 4151? For 1-1 and 1-5, the X-ray continuum (in the flare model context) could 
be dominated by a single or several strong flares located very close to the center black hole, 
while the inner disk radius from the disk line fits (7.4 and 8.3 R 9 ) gives strong supports to 
this model. For 1-2, the X-ray continuum seems to be dominated by flares far away from the 
center black hole (Rj = 17 R s ). For 1-3 and 1-4, the X-ray continuum could be dominated 
by a single or several strong flares, which are ~ r g above and coroating with the disk at 
around 14 R 9 , with only the materials in the ring around 14 R g being illuminated. The large 
dramatic decreasing in the line EWs from 1-3 to 1-4 might be due to different degree of the 
anisotropy in the X-ray continuum, i.e., larger fraction of the X-ray continuum received by 
the disk in 1-3 than in 1-4. However, it requires some fine tune in order to keep the observed 
continuum more or less constant. Alternatively, Misra (2000) suggested that the changes in 
the disk structure may explain the observed line profile variations, but, a disk undergoing 
structure changes is also hard to keep the same radiative spectrum. Finally, the decreasing 
of the fluorescent yielding may be due to the increasing ionization in the disk skin. Under 
flare context, large change in the ionization structure of disk skin is possible even for quite 
moderate variations in the observed continuum flux if the flares have different heights If 
this is case, we would expect that the flare is lower in the 1-4 than in 1-3. This should also 
reflect in the line profile. However, the uncertainty in the disk parameters do not allow us 
to verify this. 

The observed properties strongly suggest that the rapid line variations may not be due 
the variations of X-ray continuum source. Instead, it may be due to changes in geometry 
or ionization level of the reflecting medium and at the same time, these changes should not 
affect the X-ray producing medium. If so, except for the variations of the X-ray continuum 
source, what can lead to the changes in geometry or ionization level of the reflecting 
medium? 
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An absorption line at 5.8 keV is marginally detected in 1-4, but not significant in 
other intervals. Note a variable absorption line at the same energy has been discovered 
by Nandra et al. (1999) in NGC 3516. Nandra et al. (1999) interpreted this line as the 
redshifted Fe K resonant absorption line. The large redshift is either caused by the infalling 
absorbing material or by the strong gravity in the vicinity of the black hole. Given the fact 
that the line is highly variable and the large redshift, the agreement of the redshifts in this 
two objects perhaps is not purely coincidence, but have some physical origin. If so, the 
gravitational redshift explanation would be favored. 

This work is supported by Chinese National Natural Science Foundation, Chinese 
Science and Technology ministry and Foundation of Ministry of Education. 
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Fig. I. — Background-subtracted light curve (0.4 - 10.0 keV, SIS) of NGC 4151 for the long 
ASCA observation taken in May 1995. The light curve is binned using 256 s time intervals. 
The background were estimated from source-free regions. The five time intervals are also 
shown here. 
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Fig. 2. — The Fe Ka line profiles for the five time intervals, see text for details. 
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Fig. 3. — The line profile of 1-3. dotted line: best-fit two gaussians; solid line: best-fit 
disk-line model. 
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Table 1: Dual absorbed power law fits for the 1.0 - 10.0 keV continuum of the five time 
intervals 



intervals(ks) NHl NH2 cover factor index flux a % 2 /dof 



0. 


- end 


3.inii 


7 44+1-03 


0-62l°;°t 


1 Q4+0-06 


2.292 


638/570 
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- 42 


3.291°;!° 


7.92l|95 


0.60l°;l° 


i 9 q+0.22 
l- zo -0.19 


2.123 


457/441 


42. 


- 75. 


o o fi +0.49 
°- ZD -0.51 


10-17±|;i 


n 7 q+o.05 

u - ' °-0.07 


1 ™+ - 25 

l.OU_ 22 


2.322 


485/434 


75. 


- 122. 


3-06lg;tf 


6.731?;^ 


0-691°;°! 


1 qe+0.15 


2.424 


477/452 


122. 


- 168. 


q /(q+0.33 
°- 4,3 -0.37 


10.061^ 


0.491°;^ 


1 Q4+0-24 
i -' 3 ^-0.19 


2.810 


415/456 


168. 


- end 


z -^°-0.39 


6-06±i:g 


n 7 q + 0.06 

u - '^-0.06 


1 qi+0.13 
i -> 5 - l -0.13 


2.409 


411/449 



10 10 erg/cm 2 .s (flux of 2.0 - 10.0 keV, absorption corrected) 



Table 2: Two Gaussians fits for the Fe Ka line of the five time intervals. The parameters are 
central energy (E c ), width (a) and equivalent width (EW) for two Gaussians. The equivalent 
widths are calculated relative to the X-ray continuum at 6.4 keV 



intervals 


E c l (keV) 


al (keV) 


EWl (eV) 


E c 2 (keV) 
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X 2 /dof 


SIS 
I-l 
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40 it 7 7 i 
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Table 3: Disk line fits for the Fe Ka line of the five time intervals 



intervals 


Ri (Rg) 


R (Rg) 


EW(eV) 


q 


X 2 /dof 


SIS 


o o+0.8 
O-O-0.9 


1000 / 


304±^ 


-4.418:? 


1201/1119 


1-1 


7.4±8:| 


1000 / 


375lg 


-4.2l°| 
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1-2 


16.81^ 


1000 / 


1481^ 


o q+0.6 
- z -°-0.9 
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1-3 
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-2/ 
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1-4 


9.111.1 
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1711U 


-2/ 
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1-5 


s.sli:} 
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300l| 
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frozen because the fit can not give sufficient constraints on the parameters. 



